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Abstract 

Departures of thermodynamic properties of three-dimensional superfluid ^He 
from the predictions of BCS theory are analyzed. Attention is focused on devi- 
ations of the ratios A(T = {))/Tc and [C^(Tc) - Cn{Tc)] /Cn{Tc) from their BCS 
values, where A(r = 0) is the pairing gap at zero temperature, Tc is the critical 
temperature, and and C„ are the superfluid and normal specific heats. We at- 
tribute these deviations to the momentum dependence of the gap function A(p), 
which becomes well pronounced when this function has a pair of nodes lying on 
either side of the Fermi surface. We demonstrate that such a situation arises if the 
P-wave pairing interaction V(pi , P2) i evaluated at the Fermi surface, has a sign oppo- 
site to that anticipated in B CS theory. Taking account of the momentum structure 
of the gap function, we derive a closed relation between the two ratios that contains 
no adjustable parameters and agrees with the experimental data. Some important 
features of the effective pairing interaction are inferred from the analysis. 

1 Introduction 

Non-Fermi-liquid (NFL) behavior in normal states of strongly correlated Fermi systems, 
as reflected in discrepancies between experimental data and predictions of Landau's orig- 
inal quasiparticle picture P, has been the subject of intense debate during the past 
decade. To deal with superfluid states of fermionic systems, the Landau picture - involv- 
ing immortal quasiparticles and no damping - has been joined with BCS theory to form 
the standard Landau Fermi-liquid approach. This approach is adequate for describing 
conventional ("low-temperature") superconductors. However, significant deviations from 
its predictions have been observed in experiments on strongly correlated Fermi systems. 
This evidence of NFL behavior has received little attention, despite the fact that it has 
been seen at extremely low temperatures where the Landau-BCS theory should be most 
effective. 

Conspicuous examples of NFL behavior of superfluid states are found in three-dimen- 
sional (3D) liquid ^He. Specifically, in the B-state at the melting point, the recently 
measured ratio ta = A(T = 0)/Tc = 1.99 ± 0.05 [2. of the zero-temperature gap value 
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A(r = 0) to the critical temperature Tc exceeds the famihar BCS value 1.76 by 13%. 
Additionally, the ratio rc = [Cs(Tc) — Cn(Tc)] /Cn(Tc) of the normal-superfluid jump in 
specific heat, evaluated at the critical point, is significantly greater than the BCS value of 
1.43. This effect is especially prominent close to solidification: in the B-state, the excess 
reaches 30%, and it is even greater (~ 50%) in the A-state jSJ- Such departures suggest 
that the source of error lies in the conventional, oversimplified form of the BCS pairing 
interaction, rather than in a failure of the general Landau picture. It is our purpose here 
to explore the former, less radical alternative. 

It is an unfortunate feature of BCS theory in application to superfiuid ^He that the B- 
state always wins the energetic competition with the A-state, whereas experimentally the 
A-phase occupies a substantial portion of the phase diagram. To address this theoretical 
shortcoming, various authors [1100111131^1 have introduced modifications of the pairing 
interaction V ( "strong-coupling corrections" ) inherent to the presence of a pair condensate, 
noting that in BCS theory V is evaluated in the normal state. The original idea, pioneered 
by Anderson and Brinkman in Ref . |3] and later extended in Ref. |7] , has evolved into the 
"weak-coupling-plus" (WCP) theory, developed and explored by Rainer, Sauls, and Serene 
|10| El El- This theory has become a basic tool for the analysis of superfiuid phases of 
^He, including the treatment of kinetic phenomena where the quasiparticle picture 
fails. 

To describe these varied phenomena, Rainer et al. jHH [TTl [T2] have introduced about 
a dozen input parameters, composed of density- dependent weighted angular averages of 
the normal-state quasiparticle scattering amplitude. In view of the proliferation of phe- 
nomenological parameters in WCP theory, it is not well suited for the task of disentangling 
the underlying reasons for the NFL behavior of the thermodynamic properties of liquid ^He 
and the attendant failings of the standard Landau-BCS treatment. It is our contention, to 
be supported in the following analysis, that the prospects of resolving this issue are better, 
if we proceed within the more transparent quasiparticle picture. In this regard, we may 
call attention to earlier work in the same spirit in which the normal states of several 
strongly interacting Fermi systems (notably 2D liquid ^He and some heavy-fermion com- 
pounds) have been studied within an extended quasiparticle picture. The NFL behavior 
of these systems has been successfully traced to deviations of the single-particle spectrum 
from the usual Landau formula cflIp) = Pf{p — Pf)/M*, which are in turn induced by 
the momentum dependence of the effective interaction. Similarly, departures from the 
standard BCS thermodynamic relations can appear if, in the single-particle spectrum 
E{p) = \J ^^{p) + A2(j9) of the superfiuid state, the gap function A(p) acquires a strong 
momentum dependence driven by a characteristic momentum dependence of the pairing 
interaction V (cf. Refs. HH [1311211111 EE] ) . 

We note that the momentum dependence of A(p) is ignored in the WCP theory. 
Indeed, except in the arena of nucleonic pairing JH]; the momentum dependence of the 
superfiuid gap function and its empirical consequences are poorly understood. Here we 
shall attempt to rectify this situation, in the context of liquid ^He. 

In Section 2, it is shown that the momentum dependence of A(p) gives rise to devia- 
tions from the BCS thermodynamic formulas for the zero-temperature gap and specific- 
heat jump. Neglecting strong-coupling corrections, we derive a closed relation between 
these deviations that is consistent with experimental data [21 [Sj- Strong-coupling cor- 
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rections to the pairing interaction V{pi,p2), vital for reproducing the phase diagram of 
superfluid ^He 0[7|, are incorporated in Section 3. We note, however, that such mod- 
ifications of V cannot be the major factor in explaining the observed departures from 
the two BCS thermodynamic relations. For pressures P close to Pmax — 30 bar, these 
corrections are responsible for only ~ 1/3 of the excess in the specific-heat jump at the 
critical temperature, with lesser effects at lower temperatures. 

In Section 4, we return to a more detailed analysis of the momentum dependence of the 
gap function, which becomes especially strong when A(p) develops a pair of nodes located 
on either side of the Fermi surface. We demonstrate that this phenomenon occurs in the 
event that the pairing interaction V{pf,Pf) on the Fermi surface takes on a positive sign, 
rather than remaining negative as conventionally assumed. It is argued that in liquid ^He 
the change of sign is caused by the enormous enhancement of the repulsive component 
of the effective interaction with increasing density at pressures approaching the melting 
point. We find that the stronger the repulsive part of V is relative to its attractive 
component, the more tightly this pair of nodes embraces the Fermi surface, and the 
greater the anomalies in thermodynamic properties. Our conclusions are summarized in 
Section 5, and their broader implications examined. It is conjectured that the condition 
V{pf,Pf) > serves to promote D-pairing in a number of heavy-fermion systems where 
antiferromagnetic fluctuations play no significant part. 

2 Impact of the momentum dependence of the gap 
function on the BCS relations 



2.1 The B- state 

For superfluid ^He in the B-phase, non-BCS/NFL behavior manifests itself in the first 
instance through experimentally measured departures from the BCS relations 

A(0)^7r Cs{T,)-Cn{T,) _ 12 ... 

T. 7 ' a(T,) 7C(3) • ^ ^ 

for the thermodynamic ratios ta and rc appearing in the leftmost members of the two 
equations. Here ln7 = 0.577 is the Euler constant, and C,{^)i the Riemann zeta-f unction. 

We now analyze the role of the momentum dependence of the gap function A(p) in 
these deviations, neglecting strong-coupling corrections for the time being. In the B-state, 
the BCS gap equation has the conventional form 

A(p. D = - / v(p.pO *"°"^f(;;f^' A(p..r)d«. . (2) 

where dv = p'^dp/27r'^ and V{p,pi), a real function, denotes the effective P-wave pairing 



interaction. In BCS theory, the quasiparticle energy E{p) = ye'^{p) + A^(p) contains the 
single-particle energy in the normal state, given by the FL formula t{p) = pf{p — Pf)/M*, 
where M* is the effective mass. The gap function A{p) is conveniently written as a 
product 

A(p,T) = A(T)^(p), (3) 
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where A(T) = A(p = Pf,T) is the magnitude of the gap, and the factor ipip) represents 
its shape, normahzed by ipipp) = 1- In what follows we neglect the minor dependence of 
ip{p) on temperature, and assume this function obeys the equation 



Hp) = - / v(P.P.) *°°''t'f;f^'' v<PO^». 



(4) 



which is valid to a very good approximation. 

Multiplication of Eq. Q by the product 4'{p)ta.nh.{E{p)/2T)/2E{pi), integration over 
momentum p, and other simple manipulations lead to the identity 



tanh 



2Te 



tanh 



e{p) 



22] 

e2(p) + A2(T)V^2(p) 



dv = . 



(5) 



In the BCS theory developed for conventional superconductors, the pairing interaction is 
momentum-independent. Hence ip = 1, and formula (jSJ becomes 



tanh 2^ 



tanh 



BCS 



(T) 



2T 



Subtracting Eq. ^ from Eq. Q we arrive at 



dv = 0. 



(6) 



(T) 



2T 



■ (it; 



e(p) 



tanh 



v/62(p)+A2(r)^2(p) 

2T 

e2(p) + A2(T)V;2(p) 



rft; . (7) 



Departures from the BCS relations (P) will be ascribed to substantial momentum 
dependence of the gap function in the region adjacent to the Fermi surface, where the 
integration over p can be replaced by an integration over the spectrum. Thus we rewrite 
Eq. dlD as 



tanh 



BCS 



(T) 



2T 



tanh 



2T 



Ales (2^) 



e2 + A2(T) 



de 



'1 — V'^(e)) (tanh ^ — tanh ^ 



de 



+ 



r, tanh 7^ r, tanh ^ „rp — tanh ^ „rr, 
V'^(e)) ^ + V^2(e)- 2r _ 2r 



e2 + A2(T)t/'2(e) 



e2 + A2(T) 



de . 



Now let us turn to the BCS relation (1) involving A(T = 0). Setting T = in Eq. 
we have 

no r 

^ ^ t;e = /o(0) + /i(0) 



A|cs(0) 



A2(0) 



(9) 



4 



or equivalently 



where 



A2(0)-A|cs(0) 



/o(0) + /i(0) 



/o(0) 



and 



/i(0) 



A2(0) 

|e| 



■de 



A2(0) 



e2 + A2(0)^2(e) 



de. 



(10) 



(12) 



The deviation of the jump in the specific heat at T = Tc from its BCS value is 
due to the momentum dependence of the gap function as T — > Tc. One has C = TdS/dT 
with S* = — [^p In n-p + {1 — Up) ln(l — rip)], so that the jump at Tc is given by 



CsiT,)-Cnin) 



1 /9A2(T ^ Tc 

2 r 



dT 



Tp'^ip) n{p)(l — n{p)) dv . 



(13) 



where n{p) = [1 + exp(e(p)/T)]^^. As usual, the corresponding BCS formula is obtained 
from this one by setting ip = 1. Upon neglecting a small contribution from the integral 
containing the product {ip'^{e) — l)n(e)(l — n(e)), we find 



Csin) - C„(Tc) - [C,(Tc) - C„(T, 



c/JBCS 



\ dT J, 



I dT 



dT 



(14) 



To proceed further, we take T Tc in Eq. (jH}. Manipulations similar to those resulting 
in Eq. (fTUI) then lead us to the expression 



7C(3) 
47r2T2 



A2(T^Tc)-A|cs(T-^Tc 



(Tc - T) 
2Tc 



me) - 1) 



1 



cosh 



-de 



2Tc 



+ ^A\T ^ Tc 



tanh 



2Tc 



2Tc\e \ cosh 



de . 



(15) 



In deriving this result, we have employed the identity jT^] 



dx 
2^2 



tanh I 

X 



1 



2 cosh 



2 X 



703) 
47r2 



(16) 



Upon eliminating Tc — T with the aid of the BCS formula 

^2 



7C(3) 



Tc{Tc - T) 



(17) 
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Eq. (fT^ is rewritten as 



A2(T ^ T, 



A2(T ^ Tc 



/o(Tj + 



47r^ 
7C(3) 



with 



and 



oo 

— rvn ^ r. 



tanh 



2Tc 



de 



1 



2Tf.|e| cosh 



2T^ 



:i8) 



(19) 



(20) 



Now, since the Lh.s. of Eq. ()18|) coincides with the r.h.s. of Eq. ()14p. we can make the 
connection 



Cs{T,)-Cn{n) 



^ = /o(T.) + ;^A(r.). 



(21) 



It is seen that the integrands of the integrals /o(0) and /o(Tc), as well as the last term 
in the integrand of Ji(Tc), become exponentially small for \e\ ^ Tc. At the same time, 
one sees from Eqs. fjl2|) and (j^UI) that the leading terms in the integrands of /i(0) and 
/i(Tc) fall off only as e~^, at least in the interval where |'?/'(e)| < 1. Consequently, /o(0) 
and Iq{T)c receive their overwhelming contributions from the region |e| < 2Tc, implying 
that these integrals are proportional to T"^. On the other hand, both of the integrals 
/i(0) and Ii{Tc) contain an additional logarithmic factor ~ ln{Q ly/Tc) coming from the 
energy region < |e| < Qd, where fi^) is the Debye frequency. We also note that for 
|e| > A(e) = A(0)|?/'(e)|, the gap A can be neglected in the denominators of the integral 
ffT^ . Accordingly, we can write 



/i(0)/A2(0)^/i(T,)/T,^ 



l-^4(e) 



tanh —^-de 

2Tr 



(22) 



Insertion of this relation into Eq. (j2H) leads finally to 



Cs{T,) - Cn{T,) - [Cs{T,) - Cn{T,)] 



BCS 



47r2 



'A2(0) 



7C(3)A2cs(0) 



A2(0) 



(23) 



Thus, we have derived a closed relation between the departures of the two ratios r^, and 
rc from their BCS values. Importantly, beyond the thermodynamic quantities being 
connected, the relation contains no input parameters. 

The existing experimental data only allow us to test this relation for the B-state, and 
then only at P = Pmax, where A(0)/Abcs(0) — 1.15 0. Upon substituting this result 
into Eq. along with Abcs(0)/Tc = 1.76, the calculated value of the excess in the 
specific-heat jump at Tc is found to agree rather well with the experimental value of 30% 
p. 
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2.2 The A-state 



For the sake of clarity, we shall now distinguish gap functions in the A and B phases 
by corresponding subscripts. In the A-phase, one has jOj A^(p) = A'^{p)(P{n) with 
(P{n) = (3/2) sin^ 9, i.e., the gap function depends not only on the absolute value of the 
momentum p, but also on its direction n. In this case, the BCS gap equation ^ takes 
the form 

Aa(.) = - /v(.,.0^-(nO ^^"yjP^y'^^ AA(.0^, (24) 



2E(pi 



with E{p) = \j ^'^{p) + A^(p)c/2(n) and Aa(p) = ip{p)/S.x{T). The shape factor ip{j)) again 
obeys Eq. (jH), independent of the structure of the gap function. 

Repeating the same manipulations as used to reach Eq. (0), we now obtain 



tcLlltl 2^ 



(fin) 



tanh 



^e2+A^(T)d2(n),/,2(e) 
2T 



e2(p) + Ai(T)d2(n)^2 



dedn 
Air 



(25) 



In the standard BCS theory, ip = 1 and one obtains instead 



tanh ttIt 



2Tc 



tanh 



2T 

+ Ai(T)rf2(n) 



dedn 
An 



(26) 



Upon approach to the limit T Tc where A(T) vanishes, the denominator in Eq. ()26p can 
be expanded about A = 0; the integration over angles separates and is freely performed. 
Then, comparing with Eq.®, we come to what is arguably the most problematic formula 
of the BCS theory of superfluid ^He jHl El , namely 



dT 



dT 



(27) 



with K, = 6/5, leaving no room on the phase diagram for the A-phase - a prediction in 
conflict with experiment. 

Returning to Eq. and considering temperatures close to T^, we see that repair 
of the BCS phase diagram of superfluid ^He cannot be achieved by only incorporation 
of the momentum dependence of the shape factor tp{p), because of the separation of the 
integrations over the direction and over the magnitude of the momentum p. This task is 
known to be the prerogative of strong-coupling corrections j3] to the free energy of the 
superfluid state. 



3 Inclusion of strong-coupling corrections to the pair- 
ing interaction 

The origin and importance of strong-coupling corrections, proportional to A'^{T) and 
reflecting the alteration of the pairing interaction V in the superfluid state, have been 
elucidated in many articles and books jHIZllHlini- In strongly correlated Fermi systems, the 
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rigorous evaluation of these corrections, e.g. through the summation of parquet diagrams, 
is impractical, so we treat their magnitude as a phenomenological parameter. Thus, for 
T = Tc we write 



f \.Zi-k =-" + ^a(T.). (28) 

I Jc 

A salient feature of the problem is the existence of a relation between the quantities 
5a (7c) and Sb{Tc), stemming from the structure of the order parameters in the B- and 
A-phases. Correspondingly, 6a{Tc) ~ 36b{Tc), and the same ratio holds at T = |7|. 
Furthermore, the value of 5a (^c) must exceed 3(k — 1)/2 ~ 0.25 in order to outperform 
the factor n in Eq. (jTfj) and protect the A-phase from extinction. On the other hand, at 
pressures P ~ Pmax, the experimental values of the A- and B-phase specific-heat jumps 
are the same within 10%, indicating that Sb{Tc) — 0.1. Based on these considerations, 
the inclusion of strong-coupling corrections alone (without accounting for the momentum 
structure of V) is incapable of explaining the observed departures from the BCS relations 
(^. To wit: near Tc these corrections provide for only 1/3 of the measured deviations, and 
according to Ref. [7j their impact declines as T ^ 0. Thus we conclude that when one fo- 
cuses on thermodynamic properties of 3D superfluid ^He within the quasiparticle picture, 
the shape factor corrections are needed to resolve conflicts between the experimental and 
theoretical values of the ratios and tq appearing in Eq. (0). 



aA|(7e) 

dT 



dT 



1 + 5b(T, 



4 Origin of nodes in the shape factor 

Ordinarily, in studies of superfluid Fermi liquids it is tacitly supposed that the momentum 
dependence of the gap function is minor and of negligible consequence. However, the 
situation dramatically changes, if the gap function develops a pair of nodes lying on 
either side of the Fermi surface. This behavior plays a decisive role in explaining the 
departures of the ratios and from their BCS values appearing in Eq. (jT)). The 
emergence of such a pair of nodes in the shape function ip{p) does occur provided the P- 
wave pairing interaction V acquires the "wrong" (i.e., positive) sign at the Fermi surface, 
Vf = y{PF,PF) > 0. At first glance, this condition appears to rule out the existence of 
nontrivial solutions of the BCS gap equation But this is not the case. Indeed, for some 
plausible pairing interactions V used to describe strongly correlated Fermi superfiuids, it 
can happen that the sign of V{pi,P2) is positive not only on the Fermi surface, but 
everywhere in the (^1,^2) plane, yet Eq. Q nevertheless admits a nontrivial solution for 
the gap PHI- Such solutions can occur for an interaction V which, viewed in coordinate 
space, possesses a strong inner repulsive core and an outer attractive well. The strength 
of the attraction may be insufficient to compensate the repulsion in forming the pairing 
matrix elements in momentum space, but still strong enough to create a Cooper-pair 
instability. 

In what follows we employ the separation method that has been developed for efficient 
analysis and solution of gap equations [18j. The defining step of this procedure consists 
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in decomposing the interaction V(pi,p2) into a separable part and a remainder W{pi,p2) 
that vanishes when either argument is at the Fermi surface: 



V{pi,P2) = VFy^{piMP2) + W{p,,P2) ■ (29) 

The choice ip{p) = V{p,pf)/Vf meets the required condition W{pF,p) = W{p,pf) = 
for all p. Upon inserting the decomposition ()29|) into Eq. (0) one finds 



. tanh '^'^^^^ 

iP{p) = ^{p) - J W{p,p,)^^^^{p,)dv^ , (30) 

and 

r tanh ^ 

1 = -VfJ ^{p)m^^dv . (31) 

Since iP{pf) = ^{pf) = 1, the sign of the leading logarithmic part ~ ln(l/Tc) of the r.h.s. 
of Eq. (p?T|l . being negative due to Vf > 0, is opposite to that of the l.h.s. However, if the 
shape factor ip^p) changes its sign (see below), and this node occurs so close to the Fermi 
surface that the remaining part of the integral on the r.h.s. of Eq. (jHTj) becomes positive 
and overweighs unity, then nontrivial solutions of the gap equation exist • The number 
of nodes of the function ip{p) depends on the structure of the pairing potential V, and 
in particular, on the relative strength of its repulsive and attractive components, which 
governs the sign and value of the key parameter Vf- 

To illustrate and affirm these assertions, we choose a pairing interaction V(pi,p2) of 
the form ^ ^ 

V(pi,P2) = 7 ^2 I ~ 7 % I m ■ (32) 

(pi - + (pi - V2Y + /^i 

Such a model pairing potential mimics the realistic pairing interaction in superfluid ^He, 
which features a repulsive core and a long-range attractive part j2Sl 121] • To treat P- 
pairing, we evaluate the first harmonic of this interaction over the angle between pi and 
P2 and employ the result in Eq. ((201) • 

First we elucidate the emergence of the pair of the nodes in the shape factor ip{p) in 
the case of Vi? > 0. Fig. 1 shows results from solution of Eq. (jHUI) with two sets of input 
parameters: (i) (3r = Pf, Pa = O.lpF, VaN{0) = -3.0, and V;A^(0) = 45, where iV(0) = 
PfM* /tt^ is the density of states; and (ii) Pr = Pf, Pa = 0.05pF, VaN{0) = —5.5, and 
VrN{0) = 90. In both of these cases, the value of the dimensionless coupling parameter 
N{0)Vf — 2.0 is positive and quite large. 

The lower panel of Fig. 1 shows (p{p) = V{p,Pf)/V{pf,Pf), the interaction profile in 
the P-pairing channel. This function has the shape of a wide hump, with a narrow dip at 
p = Pf- The P-wave shape factor ip{p), which vanishes for p — > 0, is shown in the upper 
panel. It exhibits a narrow, sharp peak at the Fermi surface that creates two nodes, one 
on either side of the Fermi surface and lying very close to pi;'. 

Now we address the question: how does the number of nodes in the shape factor ip{p) 
change as is varied while holding Va fixed? Evidently, if the repulsion is weak, then the 
shape factor ip{p) has no nodes at all. However, if the repulsion begins to rise, such that 
Vf attains a critical value Vif, the first node of ip{p) emerges at some momentum po > pf- 
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P/Pp 

Figure 1: Shape factor ipl^p) (top panel) and interaction profile (p{p) 
(bottom panel) versus p/pr- The solid (dashed) line corresponds to the 
case Pa = 0.1 Pf, VfN{0) = 2.0 {pa = 0.05 pp, VfN{0) = 2.5). In both 
cases Tc ~ 10~^e^. The inset depicts the behavior of the shape factor 
in the vicinity of the Fermi surface. 

Whether the sign of Vif remains negative or already becomes positive will depend on the 
particular form of the pairing interaction and on the input parameters that specify it. 
With further increase of Vp, the nodal location pq moves toward the Fermi surface from 
the outside. However, the destination pq = pp turns out to be unattainable [TH] . 

Meanwhile, as Vf continues to increase and reaches a critical value V2F, a second node 
of ip{p) emerges at some momentum value below pp and moves rapidly toward the Fermi 
surface from the inside. As we shall see, the two nodes crowd the Fermi surface more and 
more tightly as V rises higher and higher above the threshold V2f- 

Fig. 2 presents results from numerical calculations of the shape factor iplp), based again 
on the the interaction V defined in Eq. ()32j). However, we now choose VaN{0) = — 1.0 and 
illustrate the sensitivity of the results to the potential parameters. The choice Pa = O.lpp 
is maintained, but Pr is increased to 2pp. From this figure, we may infer that the first 
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0.0 0.5 1.0 1.5 2.0 2.5 

Figure 2: P-pairing shape factor ijj{p) as a function oi p/pp. The inter- 
action form is assumed, with VaN{0) = 1, Pa/PF = 0.1, Pr/PF = 2, 
and five different values of the strength parameter Vr of the repulsive 
component. The curves are labelled with the corresponding values of 
the P-wave coupling parameter vp = VfN{0). 

node of iplp) emerges at some large momentum p ~ oo when vp = VifN{0) ^ —1.0. 
Also, it is important to note that in contrast to the situation found in Fig. 1, the second 
node of the shape factor now comes on the scene while Vp is still negative, with a value 
V2fN{0) ~ —0.3. Thus, we conclude that the sign of the critical value of the parameter 
Vp for emergence of the second node of ip{p) depends nontrivially on the specific choice 
of input parameters. 

In Fig. 3, we display the trajectories of the nodes of the shape factor ^(p) under 
variation of the key parameter Vp, with the same parameter setups as in Figs. 1 and 2. 
The two nodes rapidly approach the Fermi surface from opposite sides as vp = VpN{0) 
tends to toward a critical value vpc- (Once again, as seen, this critical value depend 
crucially on the specifics of the input potential parameters.) 

5 Discussion 

A complete theory of superfluid ^He must give a quantitative description of its kinetic 
phenomena as well as its thermodynamic properties. Necessarily, then, such a theory 
takes explicit account of the damping of single-particle excitations - as, for example in 
the weak-coupling-plus (WCP) theory [3 UHl HH [12] • By virtue of the scope of WCP 
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Figure 3: Trajectories of the pair of nodes that appear in the P-pairing 
shape function ip{p) of Fig. 1 (top panel) and Fig. 2 (bottom panel), 
versus the P-wave couphng parameter vp. 

theory, its phenomenological character entails the introduction of a large number of ad- 
justable parameters, which inevitably complicates the interpretation of new experimental 
results. In contrast to this approach, our concentration on the thermodynamics of super- 
fluid ^He has allowed us to retain the conceptually simpler Landau picture jT]. In essence, 
this picture stems from the basic statement, or assumption, that the ground-state energy 
and other thermodynamic quantities are functionals of the quasiparticle momentum dis- 
tribution n{p). Seemingly innocuous at first sight, the Landau assumption leads rather 
directly to the determination of this distribution function, which can be expressed in the 
same form as the Fermi-Dirac momentum distribution, n{p) = [1 +exp(e(p)/T)]~^. There 
is an important distinction, however, namely that the Landau quasiparticle spectrum 
^fl{p) = Pf{p — Pf)/M* differs from the single-particle spectrum of the ideal Fermi gas, 
since it involves an effective mass M* different from the bare mass M. In fact, this for- 
mula for epL is the most vulnerable element of Landau theory as traditionally practiced. 
The spectrum eFL(p) ceases to be meaningful close to the quantum critical point (QCP) 
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Figure 4: Critical temperature (in units of Sp = pp/2M), plotted 
against P-wave (bottom axis) and S'-wave (top axis) coupling parame- 
ters vp = VfN{0), for a model interaction of form ()32|) . The solid curve 
corresponds to P-pairing, and the dashed curve, to S-pairing. For com- 
parison we show the BCS estimate (jH^ of Tc for P-pairing (dotted 
curve) and for ^-pairing (dot-dashed curve). 

where the effective mass M* diverges. To rectify the theory, one must take into account 
contributions to e(p) from terms ~ {p — Pf)^, as established in Ref. JH]- This exten- 
sion of the Landau quasiparticle picture alters the standard Fermi-liquid formulas, with 
the implication that in the vicinity of the QCP, (so-called) non-Fermi-liquid behavior is 
predicted within the Fermi-liquid approach itself. 

A parallel situation may arise in superfluid Fermi liquids. Non-BCS behavior can be 
deduced within the standard BCS approach provided a well-pronounced momentum de- 
pendence of the gap function A{p), driven by the momentum dependence of the pairing 
interaction, is properly incorporated. In liquid ^He, the bare atom-atom interaction, as 
modelled in Ref. [5(7 by a local two-body potential function in coordinate space, con- 
tains a huge repulsive core of radius ~ 2.5 A, surrounded by an attractive, long-range 
(power-law) van der Waals component. The properties of this in- vacuum interaction are 
presumably mirrored by an in-medium pairing interaction V consisting of a local repul- 
sive component at short distance and a longer-range attractive term ■ From the 
perspective offered and supported in this paper, the existence of P-wave superfluidity 
in liquid ^He implies the presence of an attractive component Va of the effective pairing 
interaction V strong enough to outperform the repulsion at distances slightly exceeding 
the average distance tq between particles. On the other hand, at large pressures, the 
repulsive part Vr of the effective interaction increases rapidly with the density p. Indi- 
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rect evidence for the latter inference is seen in the fast growth of the sound velocity at 
pressures P ~ Pmax- 

Let us suppose that the parameter V{pf,Pf) has already turned positive at pressures 
close to Pmax, rcsultiug in the appearance of two nodes in the shape factor ip{p) of the 
gap function that hold the Fermi surface in a "vice-like grip." Furthermore, to facilitate 
the analysis, let us assume that these nodes are situated symmetrically with respect to 
the Fermi surface on either side of it, so that their locations may be specified by a single 
parameter. The plots of iplp) obtained for models based on the pairing interaction (jH^ 
indicate that this assumption is not unreasonable (see Figs. 1 and 2). It becomes more 
reasonable when long- wavelength spin fluctuations, analogous to the mediating phonons of 
conventional superconductors, dominate in the attractive part of V, since in this case the 
induced pairing interaction is symmetrical with respect to departures of the quasiparticle 
momenta from the Fermi surface. 

Numerical analysis shows that the dominant contributions to the integral /i(0) of 
Eq. (|12p come from the interval [— eo,eo]. Upon neglecting the contributions from the 
remainder of momentum space and introducing a new integration variable x = e/eo, the 
integral ()12|) takes the form 



It should be emphasized that deviations of the ratios ta and tq from their BCS values 
proportional to the ratio A^(0)/eo will depend crucially on the structure the shape factor 
ipix) in the interval < x < 1. To support this statement, we calculate the value of 
the integral ()34|) for two phenomenological shape factors: (i) ipi^) = 1 — and (ii) 
iplx) = (1 — x^)". In the latter case we take a = 1.6, a value that provides an adequate 
fit of iplx) as given in Fig. 1. Calculation of the integral in ()34|) yields K{0) = 6 in the 
first case and K{0) = 50 in the second. 

We may attempt to narrow the uncertainty in the value of the parameter eo, based on 
the observation that in the normal state at pressures P close to the melting point, the 
experimental spin susceptibility x(T) coincides with the Curie susceptibility Xc(T) = p/T 
at T > 0.4 K. This finding implies that at P ~ Pmax, the bandwidth fl = eipp) — e{p = 0) 
may be estimated as 0.2-0.3 K. We then suggest that for P ~ Pmax, the value of eo, which 
must be significantly less than Q, lies in the interval 0.05-0.1 K. Inserting this estimate 
into Eq. ()33|) . along with the two values obtained for i^(0), one obtains a range of values 
for the deviation of ta from the BCS value in the interval l%-50%. 

In Fig. 4 we present results from calculations of the critical temperature Tc for P- 
pairing and S'-pairing as a functions of the coupling parameter vp = VfN{0). Once 
again, a pairing interaction of form (j32j) is assumed to have the same set of parameters 
as in Fig. 2. Here we also take the opportunity to highlight the problematic nature of 
the standard BCS formula commonly used to estimate the critical temperature for pair 



h{0) = {A'iO)/el)K{0) 



(33) 



where 




(34) 
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condensation in fermionic systems, 

wherein A^(0) = pfM*/7t'^ is the density of states, Qd is the Debye frequency, and A is 
the couphng constant, usually identified with —Vf- Fig. 4 compares critical temperatures 
obtained with this formula (curves labelled BCS(P) and BCS(S)) with those determined 
by direct solution of the gap equation taking full account of the shape dependence of the 
gap function. Evidently, the estimate ()35p is irrelevant if Vf > 0, so meaningful com- 
parisons can only be made for vp > 0. It is seen that the formula already begins 
to fail well before Vf changes its sign. In general, then, the domain where superfluidity 
exists is wider than that where Vp is negative, a fact long appreciated within the theory 
of nucleonic superfiuids (see Refs. j2H and works cited therein). In light of this 

conclusion, the findings of Kuchenhoff and Wolfie 2,4^ on the occurrence of superconduc- 
tivity in the dilute electron gas bear re-examination, since these authors determined the 
extent of the superconducting phase from the condition Vf < 0. Referring once more to 
Fig. 4, we may also note that P-pairing wins the energetic contest with S'-pairing when 
the key parameter V^ takes on positive values. 

A similar situation may exist for unconventional Z)-pairing in heavy-fermion metals. A 
conspicuous example is CeCoIus, for which the antiferromagnetic state lies well above the 
ground paramagnetic state, implying that antiferromagnetic fluctuations are irrelevant to 
pairing in the ground state. In principle, there are several sources for attraction in the 
pairing interaction V; these include optical phonons and density fluctuations associated 
with the proximity to a critical density where electron liquid ceases to be homogeneous 
P31 . We argue that the attractive part of V changes insignificantly in switching from S- 
to D-pairing, while the repulsive part drops substantially. Thus, in certain heavy-fermion 
metals, unconventional D-pairing may win the contest with conventional S'-pairing if 
there is sufficient weakening of the repulsion in the D-wave channel relative to that in the 
5*- wave channel. 

In conclusion, we have analyzed the anomalous behavior of certain thermodynamic 
properties of superfluid phases of ^He in the framework of the Landau quasiparticle pic- 
ture. Principally, this behavior consists of departures of experimental results from two 
famous BCS relations, one quantifying the ratio of the pairing gap at zero temperature to 
the critical temperature, and the other, the ratio, at this temperature, of the difference 
between normal and superfluid specific heats to the normal specific heat. We have pro- 
posed and demonstrated that, within the quasiparticle picture, these discrepancies may 
be traced to the emergence of a pair nodes in the gap function lying close to and on 
opposite sides of the Fermi surface. We have found that a quantitative explanation of the 
anomalies requires only two phenomenological parameters, one specifying the location of 
the nodes of the gap function, and the other setting the scale of strong-coupling correc- 
tions. Finally, it should be noted that our phenomenological theory has little in common 
with the (likewise phenomenological) weak-coupling-plus theory of Refs. [7| HUl HH [12] , in 
which the momentum dependence of the gap function is completely ignored. 



(35) 
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